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Crystal structure and electronic structure of quaternary semiconductors 
Cu2ZnTiSe4 and Cu2ZnTiS4 for solar cell absorber 
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We design two new I2-II-IV-VI4 quaternary semiconductors Cu2ZnTiSe4 and Cu2ZnTiS4, and systemati- 
cally study the crystal and electronic structure by employing first-principles electronic structure calculations. 
Among the considered crystal structures, it is confirmed that the band gaps of Cu2ZnTiSe4 and Cu2ZnTiS4 
originate from the full occupied Cu Sd valence band and unoccupied Ti Sd conducting band, and kestcrite 
structure should be the ground state. Furthermore, our calculations indicate that Cu2ZnTiSe4 and Cu2ZnTiS4 
have comparable band gaps with Cu2ZnTSc4 and Cu2ZnTS4, but almost twice larger absorption coefficient 
aibj). Thus, the materials arc expected to be candidate materials for solar cell absorber. 
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I. INTRODUCTION 

Main focus of current research in solar energy and ther- 
moelectric power generation is to discover and develop 
efficient materials. An ideal thin-film solar cell absorber 
material should have a direct band gap around 1.3~1.5 
eV with abundant, inexpensive, and nontoxic elements. 
Cu-base selenide and sulfide are considered to be a 
promising class of materials for solar-cell applications^"— . 
Cu(In,Ga)Se2 (GIGS) is one of the most promising thin- 
film solar cell materials, demonstrating an efficiency 
about 20%i. However, In and Ga are expensive com- 
ponents, and the band gap is usually not optimal for 
high efficiency GIGS solar cells. Gurrently, designing 
and synthesizing novel, high-efficiency, and low cost so- 
lar cell absorbers to replace GIGS has attracted much 
attention. Among the materials that have been inves- 
tigated, the group I2-II-IV-VI4 quaternary semiconduc- 
tors Cu2ZnSnS4 (GZTS) and Cu2ZnSnSe4 (CZTSc) com- 
pounds have drawn significant interest because they con- 
tain abundant and nontoxic elements Gu, Zn, Sn, S, and 
the band gap of GZTS is about 1.5 eV— li, which is ideal 
for solar cell application. 

Motivated by the same thought, in this letter we 
report two new I2-II-IV-VI4 quaternary semiconduc- 
tors Gu2ZnTiSe4 (GZTiSe) and Gu2ZnTiS4 (GZTiS), in 
which Sn of GZTSe and GZTS is replaced by Ti. Ti is 
group IVB element and cheaper than Sn. Although its 
dioxide Ti02 has been used in dye-sensitized solar cell^, 
it has not been drawn attention to the series of group 12- 
II-IV-VI4 quaternary semiconductors yet. By employing 
first-principles electronic structure calculations, we show 
that the fundamental band gaps of GZTiSe and GZTiS 
are around 1.0'^ 1.4 eV, and the direct band gaps at F 
are round 1.2^1.6 eV, which are comparable with GZTSe 
and GZTS, thus GZTiSe and GZTiS might be useful ma- 
terials for solar cell absorber. 
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FIG. 1. (Color online) The crystal structures of CZTiSe 
(CZTiS) in (a) KS, (b) ST, (c) KS-I and (d) ST-I. Small 
yellow spheres are Se or S; middle red ones are Ti; big grey 
ones are Cu; largest black ones are Zn. 



II. COMPUTATIONAL DETAILS 



We have systematically investigated the structural and 
electronic properties of GZTiSe and GZTiS by using a 
high precise first-principles all-electron full-potential lin- 
earized augmented plane wave (FLAPW) method based 
on the density-functional theory (DFT)^. General- 
ized gradient approximation (GGA) of Perdew-Burke- 
Ernzerhof (PBE)iS is used in our calculations to obtain 
structural properties. The lattice constants a and c are 
optimized, as shown in Table I. Meanwhile, the atomic 
positions are relaxed until the force on each atom is 
smaller than ImRy/a.u. The self-consistent calculations 
are considered to be converged only when the integrated 
charge difference between input and output charge den- 
sity is less than 0.0001. We use at least 900 fc-points 
for Brillouin-zone integrations. As the GGA underesti- 
mates the band gaps, we have performed DFT-I-U cal- 
culations by using Engel and Vosko GGA plus an on- 
site Goulomb interaction Ud (EVU) to generate accurate 
band gap. Such an EVU method can correctly describe 
the band gap of d-electron zinc-blende semiconductors 
GuInSe2, ZnO, GZTS and GZTSeli^. In aU EVU cal- 
culations, the electronic correlation effects is considered 
with Ue//=4 eV for Gu, Zn and Ti, respectively. Further- 
more, It should be mentioned that the EVU calculations 
have a little influence on structural properties. 



TABLE I. Calculated lattice constant a, tetragonal distortion 
parameter rj=c/2a, GGA band gap Eg , EVU band gap 

F,f^^ and the energy difference per atom 3E relative to the 
lowest-energy structure of different chalcogenides. 



Sample 


Lattice 


a (A) 7? 


5E{meV) 


EGG-4(eV) 


Er^(eV) 


CZTiSe 


KS 


5.674 1.005 





0.56 


1.24 


CZTiSe 


ST 


5.696 L015 


2.1 


0.38 


1.06 


CZTiSe 


KST 


5.653 0.994 


34.4 


0.50 


1.17 


CZTiSe 


ST-I 


5.684 1.010 


17.3 


0.31 


0.97 


CZTiS 


KS 


5.381 0.998 





0.68 


1.42 


CZTiS 


ST 


5.418 1.018 


1.8 


0.56 


1.24 


CZTiS 


KST 


5.372 0.993 


72.4 


0.64 


1.32 


CZTiS 


ST-I 


5.403 1.011 


40.8 


0.49 


1.14 



III. RESULTS AND DISCUSSION 

The I2-II-IV-VI4 quaternary semiconductors have two 
fundamental crystal structures: one is kesterite struc- 
ture (KS) [space group 14, Fig. 1(a)], the other is stannite 
structure (ST) [space group I42m, Fig. 1(b)]. We fur- 
ther consider two structures derived from KS and ST in 
which the position of Ti and Zn are exchanged every two 
planes, named as KS-I [space group P4, fig. 1(c)] and 
ST-I [space group P42m, fig. 1(d)]. In all of structures, 
Se or S (group VI) atom is surrounded by two Cu (group 

I) atoms, one Ti (group IVB) atom and one Zn (group 

II) atom, therefore the octet rule is obeyed. 

The calculated properties of CZTiSe and CZTiS of the 
four structures are listed in Table I. We find that KS com- 
pounds have the lowest total energy for both CZTiSe and 
CZTiS, whereas ST ones have lower total energy. Other 
two structures KS-I and ST-I should be unfavorable be- 
cause of higher total energy. It seems that Ti and Zn are 
more likely to be tetrahedral coordination but not square 
coordination. Although the KS possesses a lower energy, 
the energy difference between KS and ST is small, about 
2 meV/atom. It is similar with CZTSe and CZTS that 
KS is more stable than ST only with a tiny energy dif- 
ference. Thus, it indicates that KS and ST ordering may 
coexist in real material, and the disorder effects should 
be similar with CZTS and CZTSe as discussed in Ref. 6. 

The calculated band gap of CZTiSe is smaller than 
CZTiS in the same crystal structure, which is consis- 
tent with the common expectation that the selenides have 
smaller band gaps than the corresponding sulfides. From 
the EVU calculations, estimated energy gaps Eg arc 1.42 
eV for KS CZTiS, 1.24 eV for ST CZTiS, 1.24 eV for KS 
CZTiSe, and 1.06 eV for ST CZTiSe, respectively. The 
KS ones have the largest fundamental band gaps among 
the different cation configurations. The gap values are 
found to increase more or less linearly upon increasing U 
value. For example, if Ti has Ue//=3 eV, the band gap 
wiU decrease about 0.05 eV. Thus, the calculated band 
gaps have an error bar with order of magnitude 0.1 eV. 
It should be important to compare the theoretical results 




FIG. 2. Band structure of the (a) KS CZTiS, (b) ST CZTiS, 
(c) KS CZTiSe and (d) ST CZTiSe, calculated with the 
DFT-f U method. The Fermi energy is set at zero. The high 
symmetry lines is along T (0 0) - H (1 0) - N (i | 0) - F 

(OOO)-P(HI)- 



with experimental data in the future. Nevertheless, as 
these U are typical values used in the literature^, it would 
be reasonable to estimate the gap of CZTiS and CZTiSe. 
Furthermore, the effect of different cationic occupation 
on the band gap is indicated by comparing four different 
crystal structures KS, ST, KS-I and ST-I. Within EVU, 
we find the band gap of KS-I CZTiS is 1.32 eV, whereas 
that of ST-I CZTiS is 1.14 eV. Likewise, substitution 
positions induce lattice constant to change and the lat- 
tice constant in turn affects energy band gap closely. As 
there is no available data of lattice and gaps experimen- 
tally, it is highly expected to identify structural proper- 
ties of those materials, and further control disordering 
and atomic concentration to adjust band gaps. 

The band structures of KS CZTiS, ST CZTiS, KS CZ- 
TiSe and ST CZTiSe are shown in Fig. 2, respectively. 
CZTiS and CZTiSe show rather similar band structures 
and all of them have indirect band gaps. The calculations 
indicate that the highest occupied band are mainly dxy 
state of Cu while the lowest unoccupied band are maily 
dsz^-i state of Ti. The minimum of direct band gap is 
found on the F - H line. It is 1.6 eV and 1.4 eV for KS 
CZTiS and CZTiSe at F, respectively. While the mini- 
mum is at H point for ST compounds, it is 1.3 eV and 
1.2 eV for CZTiS and CZTiSe, respectively. The different 
gap behavior should be due to different local coordina- 
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FIG. 3. The total DOS and partial DOS of KS CZTiS. 



tioii of second nearest neighbor cations: In KS lattice, 
both Cu are surrounded by four tetrahedral coordinated 
Cu, however one is surrounded by square coordinated 
Ti and tetrahedral coordinated Zn, whereas the other 
is surrounded by square coordinated Zn and tetrahedral 
coordinated Ti. While in ST lattice, every Cu has four 
square coordinated Cu, but tetrahedral coordinated Zn 
and Ti. Thus band structures have different dispersion 
relationships. 

Density of states (DOS) of KS CZTiS is shown in Fig.3. 
It is found that below -4 eV and above 4 eV, there are 
mainly S 3p states hybridized with s states of Zn, cor- 
responding to bonding and anti-bonding states, respec- 
tively. From -4 eV to -2.5 eV, there arc bonding states 
which consist of S 3p states hybridized with Sd states of 
two Cu, while from -2 eV to -0 eV, there are mainly anti- 
bonding states which consist of S 3p states hybridized 
with Bd states of two Cu. The "p-d bonding and anti- 
bonding states are separated by a narrow gap about 0.5 
eV. Beyond Fermi level, there is a peak of Ti "id states 
up to 3 cV. The calculations indicate that the band gap 
originates from the full occupied Cu d^y valence bands 
and unoccupied d'^^'^-x conducting bands. 

As shown in Table I, it indicate that KS and ST CZTiS 
(or CZTiSe) have comparable energy gaps: Eg(KS)- 
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FIG. 4. (Color online) Absorption coefficients of KS CZTiS 
(balck solid lines), ST CZTiS (red dash lines), KS CZTiSe 
(balck dot lines) and KS CZTiSe (red dot-dash line) demon- 
strating the large band-edge absorption of these materials. 
The spectra include a 20 meV Lorentzian broadening. ST 
CZiTS and CZTiSe have similar spectra as the corresponding 
KS structures. 



Eg(ST)=0.18 cV. The discrepancy of gap size induced by 
different crystal structure is similar with that of CZTS or 
CZTSe^ii. However, comparing the band gaps of lighter 
sulfides with the corresponding structure of the selenides, 
it is found that Eg (CZTiS)-Eg (CZTiSe) is about 0.2 eV. 
The values are smaller than that between CZTS and 
CZTSe, which is about 0.5 eV-^ii. The reason is that 
the two series of semiconductors have different gap char- 
acters. The direct band gaps of CZTS and CZTSe are 
derived from the Sn s state and Cu Bd states, in which 
the position of Sn s band is controlled by anti-bonding 
between Sn s state and p states of S or Se. Whereas the 
indirect gaps of CZTiS and CZTiSe consist of Bd states 
of Ti and Cu, therefore it is much influenced by the d-d 
interaction, controlled more by U term. 

In Fig. 4, it shows the absorption coefficient a(a;) of KS 
CZTiS, ST CZTiS, KS CZTiSe and ST CZTiSe, respec- 
tively, which represents the linear optical response from 
the valence bands to the lowest conducting bands. All 
four compounds have comparable absorption, although 
with different photon energy for the onset to absorp- 
tion (i.e., the band-gap energy). It is found that KS 
CZTiS has a large band-edge absorption coefficient. No- 
tably, the calculated values is almost twice larger than 
that of CZTS and CZTSo^. That imply Ti-compounds 
might have higher light transformation efficient then Sn- 
compounds. For higher energies (about 2.5 - 4.0 cV), the 
absorptions in both CZTiS and CZTiSe are decreasing 
but it is still similar with CZTS and CZTSe, although 
the conducting band characters are different. The cal- 
culated absorption coefficient is limited by considering 
direct transitions only, and at high temperatures phonon- 
assistcd transitions will be an additional contribution 
that will somewhat broaden the spectra. Therefore, it 



is worthy to extensively investigate CZTiS and CZTiSe 
to search the possible high-efficient materials, especially 
in sulfides. 



IV. SUMMARY 

Our work demonstrates that the band gaps of CZTiS 
and CZTiSe originate from the full occupied Cu 3d va- 
lence band and unoccupied Ti 3d conducting band, which 
is quite different from the well known I2-II-IV-VI4 qua- 
ternary semiconductors CZTS and CZTSe. CZTiS and 
CZTiSe might be candidate materials for solar cell ab- 
sorber since their band gaps are around 1.2 eV'--^1.4 eV, 
which is comparable with CZTS and CZTSe, although 
we should note again that the size of gap is depended on 
the choice of U. As these materials have higher absorp- 
tion coefficient theoretically, it might be convenient for 
industrial application not only because Ti is extensively 
used and cheaper than Sn. 
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